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Abstract 

The dominant long-distance contributions to the exclusive radiative decays B — > p(uS) + 7 
involve photon emission from the light quarks at large distances, which cannot be treated 
perturbatively. We point out that this emission can be described in a theoretically con- 
sistent way as the magnetic excitation of quarks in the QCD vacuum and estimate the 
corresponding parity-conserving and parity-violating amplitudes using the light-cone QCD 
sum rule approach. These are then combined with the corresponding short-distance con- 
tribution from the magnetic moment operator in the same approach, derived earlier, to 
estimate the decay rates T(B —>■ p(u) +7). The implications of this result for the determi- 
nation of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements in radiative B decays 
are worked out. 
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1 . Exclusive radiative B decays provide a valuable source of information about the 
CKM matrix. Assuming that the short- distance (SD) contribution of the magnetic moment 
operator is dominating these transitions, one derives |I| 
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where takes into account SU (3) breaking in the decay form factors and the (nominal) 
phase space differences. With this assumption and isospin invariance one also expects: 

T(B U - p+ 7 ) = 2 T{B d -> p° 7 ) = 2 T(B d -> c^ 7 ) . (2) 

These relations have been used to convert the experimental upper bound on the ratio of 
the exclusive radiative B decays 0: 

x^ty •= 034 (9o%gl °' <3) 

into a bound |V^|/|V^,| < 0.64 — 0.76, depending on the estimate of the SU(3)-breaking 
parameters in the SD-piece [I], ^|, M. While this bound is at present not competitive with 
the corresponding bound from the CKM unitarity || and from the fits of the CKM matrix 
elements || which yield 

<0.36, (4) 



V t 



one anticipates that the increased sensitivity in the radiative B decay modes projected for 
the upgraded CLEO detector and the high luminosity B factories will allow to test the 
relationships ([!]) - (0) quantitatively. 

The troublesome point, which has been raised in a number of papers [[F] [|13|], is the 
possibility of significant long-distance (LD) contributions to radiative B decays from the 
light quark intermediate states. Their amplitudes necessarily involve other CKM matrix 
elements and hence the simple factorization of the decay rates in terms of the CKM fac- 
tors involving \Vtd\ and \Vt s \ no longer holds thereby invalidating the relationships ([!]) 
and (fj). A redeeming feature, however, is that the experiments by testing ([J) can deter- 
mine in a model-independent way the extent to which the radiative transitions in question 
are dominated by the SD contributions. In general, light quark (u, c) intermediate states 
enter in the magnetic moment transitions through the corresponding Wilson coefficient 
C7 (//) (see below), and through the transitions induced by the matrix elements of the 
four-Fermion operators. In fact, the effect of the light quarks to (p) is completely 
negligible as demonstrated through the explicit calculations O, [TH], being suppressed by 



powers 0(mi/mw) 2 and 0(mi/m t ) 2 with rrii = m u or m c f\. Estimates of the contribution 
of the four-Fermion operators require a certain non-perturbative technique and have been 
worked out using quark models and the vector meson dominance (VDM) approximation. 
The purpose of this letter is to suggest an alternative technique, which treats the photon 
emission from the light quarks in a theoretically consistent and model-independent way. 

t See also [n3| for a recent discussion of this point. 
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We then combine this treatment with the light-cone QCD sum rule approach to calculate 
both the SD and LD — parity conserving and parity violating — amplitudes in the decays 
B u ,d p{&) + 7- We find that the LD contributions are negligible in the neutral 5-meson 
decays Bd — > p(u) + 7 but they may contribute up to ±20% corrections in the decay rate 
of the charged B meson — > p + 7. We work out the modified form of the relations 
(|l|) and (0), and work out the consequences of this result for the extraction of the CKM 
parameters from exclusive radiative B decays. 

2. For subsequent use, we collect the definitions used in this work. Radiative weak 
transitions at the i?-meson scale are governed by the effective Hamiltonian 



v ub v: d (c 1 o 1 + c 2 o 2 ) - v tb v; d c c 7 fS o 7 + . . . 1 , (5) 



n = % 

where we have shown only the contributions which will be important for what follows. The 
Ox, 2 are the standard four-Fermion operators 

2 = (d a V"u a )(u p V^), (6) 

p// _ y(i 4- a nd a,j3 are color indices, the Cx, C 2 are the corresponding coefficients 
(depending on the scale /i); O7 is the magnetic moment operator 

O r = ^d V (l - 75)^6 (7) 

and is the electromagnetic field strength tensor, which we take, for the photon emission 
with momentum q and polarization to be 

F llu (x)=z(e^q fl -e^q u )e^. (8) 

The coefficient includes also the effect of the four- Fermi operators O5 and 0§. For 
details and numerical values of these coefficients, see |TB|, |T7]]. We concentrate on the 
decays, B^ — > + 7 and take up the neutral B decays Bd — > p(u) + 7 at the end. The 
SD contribution to the decay rate involves the matrix element 

(fn\0 7 \B) = ^(-2i)e^{p\da^(l - j 5 )b\B(p)) , (9) 

which is parametrized in terms of two invariant form factors 

(p(k)\da^(l- l5 )b\B(p)) = [e^(q-p)-p fl (q-e^)]-2F l s (q 2 )+te tlua ^p a q^2F 2 s (q 2 ) . 

(10) 

Here p and k = p — q are, respectively, the 5-meson and p-meson momentum and is 
the polarization vector of the p meson. For the real photon emission the two form factors 
coincide, Ff = F 2 = F s . This form factor was calculated in |l] using the light-cone QCD 
sum rules, and we shall follow this approach in this paper as well. 

* Our conventions for the 75-matrix and tensor conform to those in pl| . 
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Combining the above expressions we get the SD decay amplitude 
A short = -|^^^^ W 1wPV-i[v(?'P) -p,q u ]}-2F s (q 2 = 0). (11) 

The LD contributions of the four-Fermion operators 0\, 2 involve two possibilities for 
the flavour flow, shown schematically in Fig. 1. In this paper we consider the contribution 
of the weak annihilation of valence quarks in the B meson, Fig. la,b. It is color-allowed for 
the decays of charged B mesons, and is expected to give the dominant LD contribution in 
this case. In the factorization approximation, we write 

(pj\o 2 \b) = (p|dr^|o)( 7 |«r^|5) + (frr\dr^u\o)(o\um\B) , (12) 

and make use of the definitions of the decay constants 

(0|nr^|£?> = i Pfl f B , 

<p|dl>|0> = efm p f p , (13) 

to reduce the problem at hand to the calculation of simpler form factors induced by vector 
and axial-vector currents. The two terms in fll2]) correspond in an obvious way to the 
contributions of photon emission from the loop containing the b quark, Fig. la, and from 
the loop of the light quarks, Fig. lb, respectively. The latter subprocess involves an (axial) 
vector transition with the the momentum transfer m 2 . According to our estimates, this 
contribution is much smaller than the one coming from the first term, which is dominated 
by the photon emission from a soft u-quark and which we shall consider in detail. To that 
end, we write down 

/pfrlWlfl) = -eeW"{ -i[ gfl v(q ■ p) - Plt q v ] ■ 2F 1 L (q 2 ) + e^pV • 2if(g 2 ) 

+ contact terms I , (14) 

where the two form factors describe parity-violating and parity-conserving amplitudes, re- 
spectively. Note that we have included the factor f p on the l.h.s. of flT"4] ) to make the form 
factors dimensionless. The contact terms indicated on the r.h.s. of ([14]) depend on the 
gauge chosen for the electromagnetic field and should be omitted as they are identically 
cancelled by similar gauge-dependent contact terms contributing to the second term in 
(|l2l). Physically, they correspond to contributions of photon emission from the W boson 
and arise in intermediate steps because the factorization approximation in flI3|) introduces 
charged weak currents. Adding the contribution of the operator 0\, we obtain the LD 
amplitude to the decay B u — > p + 7 in terms of the form factors F± and F 2 L , 

eG FTr _ /„ 1 



Aion 9 = -^V ub V: d (C 2 + —C^m^e^ 

x { - i[g^(q ■ V) - V„q v ] ■ 2F 1 L (q 2 ) + e^p a q p ■ 2F 2 L (q 2 )} . (15) 
We now proceed to estimate the form factors defined in fll4l). 



3. The principal contribution to the form factors is due to the photon emission from 
the u-quark which involves both contributions of small (of order 1 / m&) and large (of order 
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1/Aqcd) distances. A consistent tool to separate them is provided by the Operator Product 
Expansion, and is most easily implemented using the external field technique [Tj|. To this 
end we consider quark propagation in the background of the external electromagnetic field 
Vp, adding to the QCD Lagrangian density an extra term 5C = —ej^ n (x)V' J '(x) where 
is the electromagnetic current. The photon emission is obtained as the linear term 
in the expansion of the relevant Green functions in the external field. The corresponding 
formalism has been worked out in detail in J18, [19], EO, |22j and in this paper we only give 



a summary of the results. 

To the leading twist accuracy, that is when the quark is propagating close to the light- 
cone, all the necessary information about the magnetic photon emission is contained in a 
single (nonperturbative) matrix element of the gauge-invariant nonlocal operator ip(0)ip(x) 
at light-like separations x 2 = P0|, : 



(V , (0)cr a/3 Pexp 



i dux^gA^- e u V,j)(ux) 



ijj{x)) F = e^x{^) / du(p 7 (u)F a(3 (ux) . (16) 



The (■ ■ -)f denotes the vacuum expectation value in the external field F. Note that the 
path-ordered exponential factor includes both the gluon field and the photon field V^. 
In what follows we drop the gauge factors, assuming the Fock-Schwinger gauge x^A^ = 
x iiVij, — 0. The normalization is chosen in such a way that J^du(f)y(u) = 1, (ipip) is the 
quark condensate and the dimensionful constant x nas the physical meaning of being the 
magnetic susceptibility of the quark condensate. In the limit of a constant external magnetic 
field F a p(x) = const., the spins of quarks in the vacuum tend to get oriented along the field 
direction, with the average spin being proportional to the quark density (ipip), the applied 



field F a p, the quark charge and the magnetic susceptibility of the medium ||18|| : 



(17) 



The value of the magnetic susceptibility is large 

X(n = 1 GeV) = -4.4 GeV" 2 

and because of this the nonperturbative photon emission is numerically very important .0 

The physical meaning of the response function 1 
ticular choice of the external field plain wave 
as 



dye iqy (0\T{j'r(y)^0)a a ^(x)}\0) = W (#>[w<, - 5^] / due^^u) . (19) 



(u) becomes transparent for the par- 



In this case (|16|) can be rewritten 



Comparing this expression with the usual definition of the light-cone hadron wave functions 
p3fl it is easy to see that 7 («) defines the distribution amplitude (photon wave function) 
describing the photon dissociation into a quark-antiquark pair with the variable u being 
just the momentum fraction carried by the quark. The shape of this distribution in general 
depends on the normalization scale for the quark fields. However, an accurate analysis 



§ A larger value x = 
point 11 = 500 MeV. 



-5.7 GeV given in the first of refs.[M corresponds to a lower normalization 
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shows that unlike meson wave functions 
form already at low virtualities fi ~ 1 GeV: 



this distribution is close to its asymptotic 



7 (ii) = 6w(l — u) 



(20) 



Corrections to (|i~6l) come from higher-twist contributions which are formally suppressed 
by powers of the deviation from the light-cone x 2 . For hadron wave functions these can 
be attributed to contributions of higher Fock-components in the wave function with a 
larger number of constituents. The case of the photon is still specific, since additional 
contributions arise related to operators involving the photon field instead of the gluon fieldf]. 
These specific contributions are more important than those involving quark-antiquark-gluon 
operators (whose photon-to-vacuum matrix elements are numerically small, see p0[ [22| for 
the list of all the contributions to twist-4 accuracy and numerical estimates) and can be 
calculated exactly in terms of the quark condensate. 

The most general decomposition of the relevant matrix element to the twist-4 accuracy 
involves two new invariant functions (distribution amplitudes) g^\u) and g^\u): 



du F a p{ux 



dugf{u) 



^{u)+x 2 g^\u) 



{i/j(0)a af3 ijj(x)} F = e^(ipip) 
For one particular projection of Lorentz indices the answer is given in 



all 



[UX 



(21) 



22 



e^{i>^) foa/s gt \ duF a p(ux) x4>j(u 



x 



[1 — u) + gluons + 0(x 



which implies 



CUM_i„(2) 



2 #7 W 



-4(1 -«) 



(22) 
(23) 



To determine the functions gQ\u) and g^\u) separately, we use the operator identity 







ip(0)a a/3 ip(x) = / vdv - — ip(0)a av x v ^(vx) - (a (3) 

JO l OXp 

To prove (|24]), consider a simpler identity 



-af)r)p 







1 d - 

vdvxn- — ip(.Q)l5' l P(. vx ) ■ 



'dXr 



dv v 2 —ip(0)a a i3ijj(vx) + 2vil)(Q)a anipivx) 
L dv 



(24) 



(25) 



which can easily be checked integrating by parts. Then substitute v(d/dv) = x^(d/dx^) 
and subtract from fl25"D the first term on the r.h.s. of (^4j). By simple algebra the difference 
can be written as 



v 2 dvi>(0) XpdpG a( } - x v dpa ar} + x v d a ap v ip(vx) 



(26) 



^These additional contributions can also be thought of as contact terms, produced by operators which 
vanish on using the equations of motion. The external field technique avoids all contact terms at the cost 
of introducing additional operators. 
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where d v = d/dx v . Since for arbitrary Lorentz vectors x a ,yp, one has 

(X ■ y)a a/3 - UpOariXr) + VaVpyXr, = -tafapX^Vplh , (27) 



this expression coincides with the second term on the r.h.s. of (|24[) . 

To use fl24"|), note that the vacuum-to-photon transition matrix element of the second 
term on the r.h.s. vanishes. Substituting the general expression (|2~lj ) for the matrix elements 
in the first term, one obtains the (exact) relation 

2^ 3 f- A 9?{u)=gf\u). (28) 

Ju V 

Solving the system of equations (^), (£2^) we get: 

9?(u) = -i(l-„)(3-u), 



9?(u) = ~\{l-u)\ (29) 



which is our final result for the photon emission to twist-4 accuracy. The contributions of the 
quark- ant iquark-gluon operators can be added using (|4]) and the expressions given in |22 . 
They are numerically small and are omitted in the numerical analysis given below. Taking 
the external field in the form of the plain wave (§) we get an equivalent representation: 

dye^(0\T{j; m (y)m^(x)}\0} = 

= e^(i/)ip}[qpg atM - q a g^} / due mqx \x<j) 7 (u) + a; 2 ^ 1 ^)] + eTp(ipip){(qx)[xpg af t - x a g Pli ] 

^ 

+ x^q^Xa - q a X(3] - x 2 [q p g a ^ - q a g^]\ f du e luqx g {2) (u) . (30) 



o 



4. The expressions given above are sufficient for the description of real photon 
emission in the case that the kinematics of the particular process ensures that the quark 
propagates near the light-cone. To make a quantitative estimate of the form factors, we 
use the modification of the QCD sum rule technique, suggested in pO], |25|, |22|, [[J. The 



essence of this approach is to avoid introduction of model-dependent wave functions of the 
B meson, replacing it by a suitable interpolation operator, and using dispersion relations 
and duality to pick up the contribution of the B meson. Following this method, we consider 
the correlation function 



T»v( P ,q) = i 2 J dxe-^J dye <w <0|TOT(3/)«(0)r„6(0)6( a :)i75«(a;)}|0> 

= i[(p • q)g^v - Pt,qATi{p 2 ) + e^ai3Paq/3T 2 (p 2 ) + 0{p v ) + contact terms , (31) 

with fixed q 2 = and (p — q) 2 = m 2 . The invariant functions 7\ and T 2 can be calculated in 
QCD at large negative p 2 , in which case the above formalism to describe photon emission is 
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fully applicable: higher-twist components in the photon wave functions give rise to contri- 
butions suppressed by powers of p 2 |j] On the other hand, the discontinuity of the functions 
Ti(p 2 ), T 2 (p 2 ) at positive p 2 is saturated by contributions of meson states: 



f B ml 2F 1 L 2 (0) 
m b m 2 B — p 2 



+ ... 



(32) 



Making the usual assumption that the S-meson contribution corresponds to the spectral 
density integrated over the "region of duality", one arrives at the sum rule for the form 
factors of interest expressed in terms of integrals of photon wave functions. In this letter 
we cannot give a detailed derivation and refer the reader to [T|] for technical discussion. 
Repeating the by now standard steps, we arrive at the sum rules: 



fBm% 
mbfp 



2Ff(0)e 



-{m\-m\)lt 



U L 



4(m 2 + ut) M) 



it/ 



H 2 



+ 



3m b 



47T S 



(2u — 1) (e u — e b )u — e b In u \, (33) 



2F£(0)e- (m B 



K-">?)/« 



1 ^ e -(u/u) m i/t e{so _ m 2 /u) f ( ^ } r {u) 

u L 



4(m 2 + ut) 



u 



H 2 



[g?{u)-g?{u)\ 



3mb 
4^ 



[e u - e b )u - e b \nu L (34) 



where we have introduced a shorthand notation u 



u and neglected m 2 compared to 



ml, t, and Sq. Here t is the Borel parameter and Sq is the continuum threshold which stands 
for the cutoff in the dispersion relation restricting the region of duality for the B meson. 
In both the sum rules, the first term in the curly bracket corresponds to nonperturbative 
contributions of photon wave function, and the two last terms are perturbative contributions 
of the photon emission from u- and 6-quark, respectively. Numerically, the nonperturbative 
contribution of leading twist dominates and is the same in both the form factors. Thus, in 
our approach the parity-violating and parity-conserving LD amplitudes are close to each 
other. For completeness, we quote the corresponding sum rule from |l| for the SD form 
factor: 



2F b (0)e- (m B 



(m 2 B -ml)/t 



1 du 



U 



-(«/«K/*0( So _ m 2/ u ) J m b <t>f(u) + um p g { l\u 



mt+ut {a) 



(35) 



where 0^ 
Sec. 4 of 



a), g± (u) and g y "'{u) are the leading- twist p meson wave functions, specified in 
It is seen that the structure of sum rules for the SD and LD form factors is 



and g± (u) 



" By contact terms we indicate contributions which do not vanish after multiplication by q^. They should 
be calculated and subtracted to separate the gauge-invariant Lorentz structures. Another possibility is to 
use the Fock-Schwinger gauge for the photon field x^V^(x) = 0, V^(0) = 0, in which case the contact terms 
vanish. 
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very similar, and many of the uncertainties (such as the dependence on fg and f p ) cancel 
in their ratio. 

Before proceeding to present our numerical results, we would like to remark that the 
calculations presented here suggest that the VDM approach underestimates the LD ampli- 
tudes for the photon emission. The difference between the VMD approach and our method 
is two-fold. First, the value of the quark magnetic susceptibility in the VDM approxima- 



tion Xvdm = ~^/ m p [01 appears to be significantly below the result in (|18|). Second, the 



p-meson wave function is more "narrow" than the photon wave function, see [p3| , |20| , |T] and 
since the decay kinematics picks up contributions with almost the entire momentum carried 
by one of the constituents, the p-meson contribution in VMD is additionally suppressed. 
Thus, we expect that the VDM-type relations between B — > p7 and B — > pp amplitudes 
receive significant ( ~ 50%) corrections from contributions of the excited states. 

5. In performing the numerical analysis we conform to the values of the parameters 
given in [1]]. We note that all the three sum rules are dominated by the first term in the 
curly brackets, and since also the wave functions 7 and (p^ are similar, the ratio of the 
SD- and LD- form factors takes a simple form: 

Ft/F s ~ F 2 L /F S ~ euX{ ^ ~ 0.01 . (36) 

The only important correction to this estimate corresponds to the perturbative emission 
while nonperturbative twist-4 contributions are in fact negligible. Assuming the interval 
t ~ 5 — 10 GeV 2 for the Borel parameter, we obtain: 

Ff/F 8 = 0.0125 ±0.0010, F^/F s = 0.0155 ± 0.0010 , (37) 

where the errors correspond to the variation of the Borel parameter. Including other possible 
uncertainties, we expect an accuracy of the ratios in ([37]) to be of order 20%. This can be 
combined with the result of [I] 

Fl_ n = >/2F| w = V2Fi d ^ = 0.24 ± 0.04 , (38) 

to extract the absolute values of the LD form factors. Since the parity-conserving and 
parity-violating amplitudes turn out to be close to each other, the ratio of the LD and the 
SD contributions reduces to a number 

47r 2 m,(C 2 + C 1 /iV c ) F^_ VubV^ 
A long /A short - ff ■ ps ■ . 



Using C 2 = 1.10, Ci = -0.235, C c 7 n = -0.306 (at the scale p = 5 GeV) JT|, 0, we get 



L/S m b C e 7 ff F* 1 ; 

Since the Wilson coefficients are scale-dependent, this ratio is in general also scale-dependent 
(unless cancelled by a compensating dependence in the form factors.) Varying the scale 
p of the coefficients in the range m b /2 < p < 2m b , an additional dependence of ±10% 
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is introduced in R^J^* 1 ^ , which, however, is smaller than the error given above. To get a 
ball-park estimate of the ratio Ai ong /A s hort, we take the central values of the CKM matrix 
elements, V ud = 0.9744 ± 0.0010 §, \V U \ = (1.0 ± 0.045) ■ 10" 2 , \V cb \ = 0.041 ± 0.004 and 
\V ub /V cb \ = 0.08 ± 0.02 §, yielding, 

\A long /A short \ B ^ = \R?^\¥fM ~ io% . (41) 

\VtdVbt\ 

A quantitative analysis of the relative decay rates in terms of the CKM parameters is 
presented below. 

The analogous LD-contributions to the neutral B decays B d — > p'-f and B d — > wy are 
much smaller, a point that has also been noted in the context of the VMD and quark 
model based estimates. In our approach, the corresponding form factors for the decays 
Bd p(<^)7 are obtained from the ones for the decay B u — > fry reported above by the 
replacement e u — > e d , which gives the factor —1/2; in addition, and more importantly, the 
LD-contribution to the neutral B decays is colour-suppressed, which reflects itself through 
the replacement of the factor a± = C2 + C\/N c in (0) by a2 = C\ + CijN c . This yields for 
the ratio 

1^£_ = ^^_0.13±0.05, (42) 

where the numbers are based on using 02/01 = 0.27 ± 0.10 |2"7fl . This would then yield 

at most R B L d /^ pl ~ R^s^ = °- 05 > which in turn g ives Af^ pl /A^Cr < 0-02. We note 
that the LD-contribution just discussed is of the same order as the one expected from the 
diagrams in Fig. lb,c which we do not consider in this paper. In what follows we shall 
therefore neglect the LD-contribution to the neutral B decays. Restricting ourselves to the 
colour-allowed LD-contributions only, the relation (0) gets modified to: 



T(B U - pry) T(B U - pry) 



2T(B d -> fry) 2T(B d -> c 7 ) 



\r \r* 

^ _|_ j^Bu^py y ubV u , 



L/s v tb v td 

p(l — p) — p^ ~\~ 'rf 

2 • RL/sVnd-, y - ^ 2 + (Rl/s) 2 V^ . (43) 

where Rl/s = -Rf/s^ 7 an d Pi V are the Wolfenstein parameters p6| . The ratio T(B U — > 
p r y)/2T(B d — > p 7 )(= T(B U — > p r y)/2T(B d — ► w 7 )) is shown in Fig. 2 as a function of 
the parameter p, with rj = 0.2, 0.3 and 0.4. This suggests that a measurement of this 
ratio would constrain the Wolfenstein parameters (p,rj), with the dependence on p more 
marked. In particular, a negative value of p leads to a constructive interference in B u — > p 7 
decays, while large positive values of p give a destructive interference. This behaviour is in 
qualitative agreement with what has been also pointed out in |5[]. 

The ratio of the CKM-suppressed and CKM-allowed decay rates in (|l|) for charged B 
mesons likewise gets modified due to the LD contributions: 

T(B U -> fry) 2 2 2 

T(B -> i^* 7 ) = +r/] 

x {l + 2 " ^K ^j 1 ,"^"^ + (^/5) 2 K 2 d (1 f , (44) 



9 



where A = 0.2205. Using the central value from the estimate £ u = 0.59 ± 0.08 JTJ we 
show in Fig. 3a the ratio ([44]) as a function of p for 77 = 0.2, 0.3 and 0.4. It is seen that 
the dependence of this ratio is rather weak on the parameter r] but it depends on p rather 
sensitively. To show the effect of the LD contribution, we compare in Fig. 3b the predictions 
for the ratio (144]) as a function of p with and without the LD-contribution, fixing rj = 0.3, 
the central value for this parameter obtained from the CKM fits ||. It is seen that the 
effect of the LD contributions on this ratio is small, and is comparable to ~ 15% uncertainty 
in the normalization due to the S'f/(3)-breaking effects in the form factors. Neglecting the 
colour-suppressed LD contributions, as argued above, the ratio of the decay rates for neutral 
B meson is not effected and to a good approximation the SD-dominated result 

T nBZTj =^-<>?W\ (45) 

still holds. Finally, combining the estimates presented here and in Jl| for the form factors 
and restricting the Wolfenstein parameters in the range —0.4 < p < 0.4 and 0.2 < 77 < 0.4, 
as suggested by the CKM-fits ||, we give the following range for the absolute branching 
ratios: 

B(B U -> pry) = (1.9 ± 1.6) x 10~ 6 , 

B(B d -> pj) ~ B[B d -> CJ7) = (0.85 ± 0.65) x 10~ 6 , (46) 

where we have used the experimental value for the branching ratio BiB — > K* + 7) = 
(4.5 ± 1.5 ± 0.9) x 10~ 5 [pgfl , adding the errors in quadrature. The large error reflects 
the poor knowledge of the CKM matrix elements and hence experimental determination 
of these branching ratios will put rather stringent constraints on the parameters p and 
77. Similar constraints would also follow from the measurement of the inclusive radiative 
decays B — > X d + >y |2J|, for which a branching ratio B(B — > X d + >y) = (1.0 ±0.8) x 10~ 5 is 



estimated using the measured inclusive rate BiB — > X s + 7) = (2.32 ± 0.57 ± 0.35) x 10 



-4 



To summarize, we have presented a QCD sum rule-based calculation of the contribution 
of the weak annihilation to exclusive radiative rare B decays, which is expected to dominate 
the LD contribution for the decays of charged B mesons, B^ — > p^ + 7. The numerical 
effect of the weak annihilation amplitude depends on the value of the CKM- Wolfenstein 
parameter p, and is typically estimated in the 10% range for the presently allowed values of 
the parameters p and 77. The corresponding LD contributions in the neutral 5-meson decays 
B d — > p(w)7 are found to be insignificant. The interference of LD and SD contributions 
provides a possibility to measure the sign of the Wolfenstein parameter p from the ratio 
(|43|). The branching ratios B{B U — > /ry) and B{B d — > p{oS)i) likewise are sensitive to the 
sign and magnitude of p, increasing with large negative values of p. Our investigations 
strengthen the expectations that exclusive radiative rare B decays, very much like their 
inclusive counterparts, are dominated by SD contributions, which in the context of the 
standard model implies that these decays are valuable in quantifying the CKM parameters. 
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Figure Captions 

Figure 1: Weak annihilation contributions in B u — > p7 involving the operators 0\ and O2 
denoted by ® with the photon emission from a) the loop containing the b quark, b) the 
loop containing the light quark, and c) the tadpole which contributes only with additional 
gluonic corrections. 

Figure 2: Ratio of the neutral and charged 5-decay rates T(B U — > p r y)/2T(B d — > pj) as a 
function of the Wolfenstein parameter p, with r\ = 0.2 (short-dashed curve), 77 = 0.3 (solid 
curve), and r\ = 0.4 (long-dashed curve). 

Figure 3: Ratio of the CKM-suppressed and CKM-allowed radiative 5-decay rates 
T(B U — > p r y)/T{B — > K*ry) (with B = B u or Bj) as a function of the Wolfenstein parameter 
p, a) with i] = 0.2 (short-dashed curve), r\ = 0.3 (solid curve), and i] = 0.4 (long-dashed 
curve); b) comparison of the result obtained by neglecting the LD contribution (dashed 
curve) with the one including the LD contribution (solid curve), both evaluated with r\ = 
0.3. 
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